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Measurement of the transmitted intensity from a coherent monomode light source through a series 
of subwavelength slit arrays in Ag films, with varying array pitch and number of slits, demonstrate 
enhancement (suppression) by as much as a factor of 6 (9) when normalized to that of an isolated 
slit. Pronounced minima in the transmitted intensity were observed at array pitches corresponding 
to Aspp, 2 Aspp and 3Aspp where Aspp is the wavelength of the surface plasmon polariton (SPP). 
Increasing the number of slits to more than four does not increase appreciably the per-slit trans- 
mission intensity. These results are consistent with a model for interference between SPPs and the 
incident wave that fits well the measured transmitted intensity profile. 
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Since the first experimental report of "extraordinary 
optical transmission" through subwavelength hole ar- 
rays ll| , considerable effort has been devoted to the essen- 
tial physics of the process in both hole [2, 3, 5, 6, 7J and 



slit |8|, [9|, ll0|, lUl, ll3| arrays. Some studies [2|, 13, l8[ have 
interpreted the transmission spectrum as excitation of 
delocalized surface plasmon Bloch modes and identified 
peaks in the transmission with wavelengths equal to in- 
teger multiples of the array pitch. Others [H BBQ 



concluded that this condition should be associated with 
transmission minima. Still other theoretical work 0, 0, @] 
interpreted the transmission line shapes as "Fano pro- 
files," involving interferences between surface and inci- 
dent propagating modes. Further spectral transmission 
measurements [15] revealed that, normalized to a single 
aperture, suppression as well as enhancement was a char- 
acteristic prop erty of hole and slit arrays, and interfero- 
metric studies [l6|, [13, EH showed that the contribution 
of transient diffracted surface modes are as important 
as the surface plasmon polariton (SPP) guided mode in 
the immediate vicinity of the subwavelength object. The 
experimental setup reported in Refs. [l|, [l5[ consisted of 
an incoherent, broad-band light source passed through a 
scanning spectrophotometer and focused on fixed-period 
subwavelength hole and slit arrays. Transmitted inten- 
sity was detected in the far field as a function of the 
scanned wavelength. In that work the spectral resolu- 



tion and coherence length of light incident on the ar- 
rays therefore depended on instrumental parameters, and 
these in turn can affect the position and shape of the 
measured spectral features. Furthermore, the frequency- 
dependence of the dielectric constant of Ag and other 
real metals is non-negligible in the range of typical wave- 
length scans from 450 to 900 nm. 

In order to remove these measurement ambiguities 
and experimentally test the various theoretical interpre- 
tations, we have undertaken a series of high-resolution 
measurements of the transmission through slit arrays in 
which the spectral source is coherent, monomode and at 
fixed frequency. The transmission measurement setup 
consists of a A = 514.5 nm, 5 mW, TEMoo light beam 
from an Ar ion laser and aligned to the optical axis of an 
inverted microscope. The beam is focused at normal in- 
cidence onto the sample surface through the microscope 
condenser and polarized TM (magnetic H-field compo- 
nent parallel to the long axis of the slits). Light intensity 
transmitted through each slit array is then gathered by 
a 50X microscope objective with a numerical aperture 
(NA) of 0.45 and detected with a liquid-nitrogen-cooled, 
charged-coupled device (CCD) array detector. Light in- 
tensity transmitted through each slit array is obtained by 
integrating the signal over the entire region of interest in 
the CCD image and subtracting the background originat- 
ing from electronic noise. Per-slit transmission intensities 
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FIG. 1: Two typical elements in the overall structure layout. 
Panel (a) shows N — 4, p — 450 nm. Panel (b) shows N = 4, 
p = 700 nm. Each slit is FIB milled through a 200 nm thick 
silver layer. Dimensions of each slit are 50 nm wide and 10 /xm 
long. 
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are obtained by correcting the transmitted intensity for 
the calculated collection efficiency of the microscope ob- 
jective lens and normalizing the transmitted intensity for 
each series of gratings to the intensity collected from a 
single-slit structure. The series of slit arrays were milled 
with a focused-ion-beam (Ga + ions, 30 keV) in a 200 nm 
thick layer of silver evaporated onto a flat fused-silica 
microscope slide. The layout of slit-array structures con- 
sisted of a matrix of 9 rows and 140 columns. Each row 
was indexed by the slit number N in the array and varied 
from N = 1 — 9. Each column was indexed by the array 
pitch p starting from the first column at p — 150 nm and 
incremented by 5 nm with each successive column. Thus 
the pitch varied over a range from p = 150 — 845 nm; 
from less than Aspp to greater than 2A$pp. Each slit 
was milled 50 nm wide, 200 nm deep, and 10 /im long 
as shown in Fig.[TJ The structured silver layer was cov- 
ered by a second microscope slide; optically contacted 
to the silver surface by index- matching fluid (n = 1.46) 
so that the index change at the dielectric-silver interface 
was identical at both the input (incident) and output 
(transmitted) planes. The transmitted intensity of each 
successive array along a given row was recorded in the 
far-field by the CCD as the sample was stepped using 
an X-Y translation stage. The results are summarized in 
Fig.El 

Taking into account the collection efficiency of the 
microscope objective and the far-field angular distri- 
bution of the slit grating diffraction modes, we define 
r] = |i^jv| 2 /|i^i| 2 as the ratio of the magnetic field in- 
tensity at the output aperture of each slit in an array 
of N slits to the magnetic field intensity at the output 
aperture of an isolated slit. Figure [2] plots r] vs. array 
pitch for all arrays N = 1 — 9 (except N = 5, omitted 
due to defective fabrication). The results show that the 
transmission intensity for all arrays exhibits very simi- 



FIG. 2: (color online) Normalized transmission intensity r} 
vs. grating pitch (in micrometers on the lower abscissa and 
normalized to Aspp on the upper abscissa) for a series of slit 
arrays N = 1 — 9. Gratings with N = 5 were omitted due 
to defective fabrication. The wavelengths Ao, Ai, and Aspp 
are, respectively, the free-space wavelength, the wavelength 
in fused silica (n = 1.46) and the wavelength of the surface 
plasmon polariton. 



lar behavior with transmission dropping to a minimum 
0.1 r] at an array pitch equal to Aspp, then rising to a 
broad maximum ~ 6 rj before repeating similar behavior 
around 2 Aspp. The position of the minima are in accord 
with earlier predictions [ll|, [l2jl and simulations [l3|, 
and at variance with others [2|, 13] predicting transmission 
maxima at Aspp. 

The wavelength Aspp was calculated from the 
usual formula for the guided surface wave [lo[ , 



nspp 
Aspp 



e„ 
Ao 



(1) 



nspp 



where Ao is the free-space incident wavelength, e m , 
are the dielectric constants of the metal and adjacent 
dielectric respectively, and nspp is the effective surface 
index of refraction. In the present experiments the 
dielectric constant of the structured silver sample was 
measured directly by ellipsometry at Ao = 514.5 nm 
and determined to be e m = — 9.3 + 0.18i. The dielectric 
constant of the fused silica substrate is = +2.13, and 
therefore A SPP = 309.5 ±0.1. 

Figure [3] plots the maximum and minimum values of 
r\ for each of the N grating series. Enhancement above 
single-slit transmission up to a factor of ~ 6 is observed 
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FIG. 3: (color online) Normalized per-slit transmission inten- 
sity maxima and minima vs. the number of slits elements in 
array N. The labels 1-3 refer to the transmission intensity 
maxima and minima from left to right shown in Fig. [2] 



as TV increases up to N = 4. At and above N = 4, 
adding additional grating elements to the array does not 
significantly enhance the transmission. Similar behavior 
is observed for the transmission minima. These measure- 
ments support the view that transmission enhancement 
is dominated by nearest-neighbor slit scattering and in- 
terference with the strength of the local interaction effec- 
tively screening contributions from more distant array el- 
ements. If Bloch surface modes, delocalized over the full 
extent of the array, played a dominant role, one would 
expect the per-slit intensity to increase with the number 
of elements in the array. 

Since the positions of transmission suppression and en- 
hancement as a function of period are essentially indepen- 
dent of N, we can analyze the mechanism responsible for 
modulation by concentrating on the simplest case N = 2. 
The normalized per-slit transmission intensity r] of an ar- 
ray of slit pairs with varying pitch p is shown in Fig. 
The intensity 77 is plotted on a linear scale as a func- 
tion of p for devices milled into a Ag film of thickness 
t = 300 nm (hollow circles) in addition to the Ag film 
of thickness t = 200 nm described earlier (solid circles, 
replotted from Fig. [3]). Comparison of the two data sets 
shows that the rj(p) modulation is essentially invariant 
with t and is therefore governed by the interaction be- 
tween the two slits mediated by surface waves running 
along both facets of the structured metal film. Periodic 
minima are measured at slit-slit distances corresponding 
to integer multiples of Aspp, (p = nXspp,n = 1,2, . . .). 
This observation is consistent with recent theoretical pre- 
dictions for a two-slit system [l4| (albeit for a structure 
with only one SPP-sustaining surface). 

We have developed a simple, first-order model for rj(p). 



Figure [4^ shows the essential idea. Two surface waves, 
one at each slit, are generated by diffractive scattering 
of a normally incident plane wave. The respective SPPs 
counterpropagate along the surface, creating a standing 
surface wave and interfering with the incident wave at the 
opposite slit. The intensity transmitted to the other side 
of the metal film is proportional to the resulting modu- 
lated intensity at each slit opening. An identical process 
takes place on the exit side of the film. Counterpropagat- 
ing surface waves are launched by diffractive scattering at 
the slit exits and interfere with the directly propagating 
mode at the opposite slit location. In Fig.3^ (3 designates 
the amplitude launching efficiency of the SPP by a given 
slit at either entrance or exit aperture. More specifically 
(3 is defined by the ratio of the launched SPP magnetic 
H-field complex amplitude to that of the incident wave 
at the slit entrance or to that of the propagated mode 
at the slit exit. Similarly (3 f designates the conversion 
efficiency of the SPP complex amplitude back to prop- 
agating modes at either side of the slit opening. This 
first-order interference process yields a net transmission 
intensity (normalized to that of an isolated slit) given by, 
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( A SU rf ) P + ^ _ } 



(2) 

In Eq.[2] A sur f is the surface wavelength, (3q = \(3\, (3' = 
\/3'\ are the magnitudes of the launch and conversion ef- 
ficiencies at the slits. The phase cp = arg (/?/?') is the 
phase associated with the SPP <-> propagating wave con- 
version, exclusive of the phase accumulated along the sur- 
face, (27r/A sur f)p. Refining the model by taking multiple 
reflections into account at the slits results in the following 
closed-form expression, 



r/ (oo) (p)= l + (/3o/3o) 2 -2/?o/?oCos 



. ( Asurf ) P ^ _ } 



(3) 

A fit of 77°° (p) to the combined set of experimental trans- 
mission data for both t = 300 nm and t = 200 nm is 
shown in Fig. [4)3 (solid red curve). Excellent agreement 
is obtained using fitting parameters A sur f = 307 nm, 
P0P0 = 0.24, and tp = tt. The best-fit value for A sur f is 
slightly less than the theoretical value Aspp = 309.5 nm. 
Considering that most of the surface wave propagation 
takes place in the "near zone" 16|, 18, 2l[ where transient 
modes contribute to form a composite surface wave, this 
result is not surprising. The shape of rj°°(p) is remi- 
niscent of the transmission characteristics of a "lossy" 
two-mirror Fabry-Perot resonator of free spectral range 
A A = A sur f, full- width at half-maximum d A = 0.29A sur f 
and finesse F — AX/dX = 3.4. The Fabry-Perot profile 
suggests the influence of multiple surface wave reflections 
at the slit sites. A plot of the first-order model rf^{p) 
is also included in Fig. [lb (dashed blue curve), using the 
fitting parameters above. The essential profile of the nor- 
malized transmission as a function of p is already well 
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FIG. 4: (color online) Panel (a) shows a schematic of the 
interference model used to fit the two-slit transmission in- 
tensity. Panel (b) shows the two-slit transmission intensity 
normalized to single-slit transmission as function of slit-slit 
separation and plotted on a linear scale. The data (filled cir- 
cles) show the transmission profile for a 200 nm thick Ag film; 
open circles show similar data for a 300 nm thick Ag film. The 
dashed line shows a fit using the first-order interference model 
of Eq.[2] (blue curve) and the solid line shows a fit using the 
infinite-order model of Eq.[3] (red curve). 



reproduced by rj^{p). The first-order fit suggests that 
the formation of transmission minima is predominantly 
controlled by interference at the slit openings rather than 
by the presence of higher order multiple reflections. The 
positioning of the minima at p = nAspp, (n = 1,2,...) 
is due to (p = 7r and is in agreement with the findings of 
Ref. [14J although the search for a simple but convincing 
physical explanation for this phase continues. 

In summary we have measured the transmitted far- 
field intensity through a series of subwavelength slit ar- 
rays as a function of array pitch and have determined that 
the minimum per-slit transmission at the array output 



facet occurs for an array pitch equal to the wavelength of 
the surface plasmon polariton. We have also determined 
that the per-slit transmitted intensity does not increase 
appreciably above an array size greater than TV = 4. 
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